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13.  ABSTRACT 


The  objective  of  this  research  was  to  provide  experimental  values  for  cone 
drag  coefficients  in  the  near  free  molecule  and  transitional  flow  regimes  including 
an  examination  of  the  effects  of  Mach  number,  wall  speed  ratio,  bluntness,  and 
specific  heat  ratio.  A  modulated  force  probe  with  a  sensitivity  in  the  microgram 
range  was  employed  to  measure  the  drag  forces  on  9°  half-angle  cones  in  the  range 
6<Mw<lQ,  l<Re  ^<50  with  an  RMS  dispersion  of  +2%.  A  knowledge  of  the  flow  properties 
in  the  free  jef.  employed  for  this  study  permitted  the  calculation  of  cone  drag 
coefficients.  The  resulting  values  cover  the  flow  regimes  from  near-free-molecule 
to  transitional. 

The  experimental  data  are  presented  as  a  function  of  the  dimensionless  parameters: 
wall  Re  number,  and  free-stream  Knudsen  number.  In  addition  an  "effective"  Knudsen 
number,  based  on  reflected  incident  molecular  collisions  employing  velocity  dependent 
cross  sections  is  shown,  based  on  limited  data,  to  provide  a  valid  correlation  for 
various  values  of  M  number  and  specific  heat  ratio  as  long  as  cone  angle  and  wall  speed 
ratio  are  fixed.  The  relationship  between  Kn  .  and  the  parameters  commonly  employed 

A  X  y  U 

in  the  literature  is  examined  and  shown  to  be  consistent  in  limiting  cases. 
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ABSTRACT 


The  objective  of  this  research  was  to  provide  experimental  values 
for  cone  drag  coefficients  in  the  near  free  molecule  and  transitional 
flow  regimes  including  an  examination  of  the  effects  of  Mach  number, 
wall  speed  ratio,  bluntness,  and  specific  heat  ratio. 

A  modulated  force  probe  with  a  sensitivity  in  the  microgram 
range  was  employed  to  measure  the  drag  forces  on  9°  half-angle  cones 
in  the  range  6  <  <  10,  1  <  Rew  ^  <  50  with  an  R.M.S.  dispersion 

of  +2%.  A  knowledge  of  the  flow  properties  in  the  free  jet  employed 
for  this  study  permitted  the  calculation  of  cone  drag  coefficients. 

The  resulting  values  cover  the  flow  regimes  from  near-free-molecule 
to  transitional. 

The  ^experimental  data  are  presented  as  a  function  of  the  dimensionless 

parameters:  wall  Reynolds  number,  Re  and  free-stream  Knudsen  number, 

w ,  x. , 

Kn  In  addition,  an  "effective"  Knudsen  number,  Kn  .  _  based  on 

09 ,  D  ’  ri,D, 

reflected-incident  molecular  collisions  employing  velocity  dependent 
cross  section  is  shown,  based  on  limited  data,  to  provide  a  valid  correlation 
for  various  values  of  Mach  number  and  specific  heat  ratio  as  long 
as  cone  angle  and  wall  speed  ratio  are  fixed.  The  relationship  between 
Knfi  D  and  the  parameters  commonly  employed  in  the  literature  is  examined 
and  shown  to  be  consistent  in  limiting  cases. 
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1.  INTRODUCTION 


Although  considerable  progress  has  been  witnessed 
during  the  last  few  years  in  providing  experimental  data  and 
theoretical  predictions  on  the  drag  of  slender  cones  at  zero 
angle  of  attack  in  the  rarefied  hypersonic  flow  regime,  a 
completely  satisfactory  understanding  of  the  problem  seems 
still  not  to  be  in  hand.  An  excellent  critique  of  near  free 
molecule  theories  and  a  discussion  of  available  experimental 
results  have  recently  been  presented  by  Keel  and  Willis  (Ref,  l) , 
It  contains  numerous  references  to  earlier  studies  which  will 
not  be  repeated  here. 

The  objective  of  this  study  was  to  provide  further 
experimental  information  on  9°  half-angle  cones  in  the  Mach 
number  range  from  6  to  10  and  to  assess  the  effects  of  specific 
heat  ratio,  bluntness,  and  wall-to-stagnat ion-temperature  ratio 
(wall  speed  ratio),  A  correlation  parameter,  the  Knudsen  number 
characterizing  reflectea-incident  collisions,  incorporating 
temperature-dependent  collision  cross  sections,  is  introduced 
and  its  usefulness  5.3  examined  by  considering  selected  sets  of 
experimental  data  from  the  literature. 
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2,  DESCRIPTION  OF  MEASUREMENT  TECHNIQUE 

This  section  describes  the  operating  principle  of 
the  modulated  force  probe  and  the  calibration  procedure 
employed  for  these  measurements. 

2  * 1  Principle  of  Operation 

Low  density  wind  tunnels  seldom  are  capable  of  pro¬ 
ducing  reasonably  large-scale  flow  fields  with  mean-free-paths 
greater  than  5  to  10  mm.  Since  it  may  easily  be  shown  that  the 
drag  force  on  a  model  with  dimensions  of  order  millimeters  in 
a  hypersonic  cold  flow  at  Kn  n  a.  0(10)  is  the  order  of  milli- 

<»  |  V 

grams,  it  is  clear  that  a  measurement  technique  suitable  for 

the  transitional  (0.3  <  Kn  _.  <  3)  and  near  free  molecule 

» ,  u 

(3  <  Knw  D  <  10)  regimes  must  exhibit  a  sensitivity  of  10  to 
100  ygrams  to  provide  a  "l%"  measurement. 

For  this  purpose,  an  instrument  was  developed  at  the 
VKI  whose  operation  is  based  on  the  principle  of  flow  modulation. 
Operational  details  may  be  found  in  Reference  2,  Briefly,  modu¬ 
lating  (turning  "on"  and  "off"  periodically)  the  flow  acting  on 
a  body  attached  to  a  typical  ceramic  force  transducer  produces 
an  A.C,  output  whose  peak-to-peak  amplitude  is  directly  propor¬ 
tional  to  the  applied  force,  A  P.A.R,  Model  JB-6  lock-in  ampli¬ 
fier  (a  highly  tuned  amplifier  -  synchronous  demodulator), 
employing  the  modulating  frequency  as  a  reference,  allows  the 
signal  of  interest  to  be  extracted  even  when  very  low  values  of 
the  s ignal-to-ncis e  ratio  exist  ,  A  schematic  of  the  instrument 
chain  is  shown  in  Figure  1, 

The  flow  field  in  which  these  drag  measurements  were 
performed  is  a  free-jet  expansion  which  provides  a  high  Mach 
number,  low  density  flow  (described  in  detail  in  the  following 
chapter)  without  the  attendent  problem  of  boundary  layer  growth 
experienced  when  the  more  common  de  Laval  geometry  is  employed. 

In  addition,  the  small  size  of  the  jet  orifice  allows  the  -flow 
to  be  modulated  quite  easily.  A  disc  with  a  semi-circular  window 
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is  positioned  very  close  to  the  orifice  plane  with  the  result 
that  the  flow  is  periodically  turned  on  and  off  at  its  source* 
A  resolver  on  the  disc's  shaft  then  provides  the  reference 
signal  for  the  lock-in  amplifier.  The  modulation  frequency  is 
15  Hz. 


The  probe  itself,  shown  in  Figure  2,  employs  a  iead- 
zirconate-titanate  ferroelectric  transducer  (PZT-5  manufactured 
by  Brush-Clevite )  to  which  the  sting  and  model  are  rigidly 
connected.  Much  of  the  sting  is  shielded  to  reduce  the  ever¬ 
present  tare  drag. 

2 . 2  Calibration  of  the  Force  Probe 

The  modulated  force  probe  is  calibrated  by  replacing 
the  model  with  a  wire  loop  and  then  inserting  a  straight  well- 
defined  segment  of  this  loop  into  the  uniform  and  known  field 
of  a  permanent  magnet.  A  precisely  measured  A.C.  currant  through 

the  wire  loop  thus  produces  an  accurately  known  axial  force 

•  , 
acting  on  the  transducer  due  to  the  effect  ,A  schematic  of 

the  calibration  system  appears  in  Figure  3® 
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3.  FLOW  FIELD  PROPERTIES 

In  this  section,  the  "ideal"  free  jet  flow  properties 
are  described,  to  present  a  simple  picture  of  the  flow  field 
employed  for  the  force  measurements.  Then  the  experimental  ar¬ 
rangement  is  discussed  with  particular  attention  to  design 
problems  encountered  when  heated  flows  are  employed.  Finally, 
the  real  or  non-ideal  free  jet  is  discussed  with  the  e.mpnasis 
on  problems  peculiar  to  this  investigation. 

3.1  The  Ideal  Free  Jet  Flow  Field 

The  free-jet  expansion  (a  highly  underexp ended  jet) 
is  employed  to  produce  high  Mach  number,  low  density  flows.  A 
schematic  diagram  of  the  jet  structure  is  shown  in  the  upper 
half  of  Figure  5.  As  it  is  a  nearly  spherical  expansion,  flow 
gradients  are  important;  however,  the  use  of  sufficiently 
small  models  circumvents  to  a  large  degree  this  problem, 

l 

The  case  of  a  jet  exhausting  isentropically  into  a 
vacuum  has  been  studied  in  Reference  ?  and  the  results  quoted 
are  now  widely  used  to  determine  "ideal"  local  flow  properties 
in  a  free-jet  flow.  The  Mach  number  or.  the  jet  centerline  may 
be  expressed  uniquely  in  terms  of  x/D‘  ,  the  distance  from  the 
orifice  to  the  point  in  question  divided  by  the  effective  nozzle 
diameter  for  one-dimensional  flow,  and  y,  the  specific  heat 
ratio.  The  Mach  number  and  density  ratio  are  plotted  versus 
x/D  in  the  love*-  part  of  Figure  5.  Thus,  the  Reynolds  number. 

Re  ,  (or  Knudsen  number  Kn  )  at  a  fixed  x/D*  can  be  altered 
independently  of  M#  simply  by  changing  the  stagnation  pressure, 
Pq.  The  test  procedure  in  a  free  jet  is  therefore  to  select  M  , 
which  fixes  x/D  ;  then  perform  the  desired  measurements  at  this 
value  of  x/D  by  varying  po  which  directly  varies  Re^  or  Kn^, 
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3.2  Experimental  Arrangement 

Two  distinct  arrangements  have  been  employed  to  pro¬ 
duce  the  desired  flow  field  depending  on  whether  or  not  the 
flow  is  to  be  heated. 

Cold  flews  (room  temperature)  emanate  from  a  50  liter 
settling  chamber  through  a  sharp-edged  nominal  5  ram  diameter 
orifice.  The  30  cm  diameter  chopper  disc  is  positioned  such 
that  the  plane  of  the  disc  is  vitnin  0.5  mm  from  the  plane  of 
the  orifice.  Stagnation  pressures  are  measured  with  a  Texas 
Instrument  Bourdon  Tube  Gauge  with  an  absolute  accuracy  of  *10 
microns  Kg,  Stagnation  temperatures  are  determined  with  thermo¬ 
couples  and  mass  flows  are  measured  with  precision  "ROTA" 
rotameters . 


Heated  flows  are  prodcued  with  an  inconel  heater  shown 
in  Figure  6,  Two  inconel  coils,  each  carrying  the  gas  load, 
bring  d.c.  current  to  and  from  the  one-liter  inconel  settling 
chamber  which  is  surrounded  by  an  auxiliary  oven.  Radiation 
shielding  protects  the  model/ transdu cer  system  from  excessive 
radiative  heat  transfer.  Measurements  of  the  stagnation  temper¬ 
ature,  Tg,  with  thermocouples  and  calculation  of  Tp  using  the 
sonic  flow  technique  -  requiring  a  knowledge  of  the  effective 
orifice  diameter  D  ,  stagnation  pressure,  pp,  and  mass  flow  m  - 
were  in  agreement  over  a  vide  range  of  operating  conditions. 
Temperatures  up  to  1000°K  car.  easily  be  produced  and  higher 
temperatures  are  certainly  possible  if  and  when  desired,  although 
the  radiant  heat  transfer  losses  then  become  significant. 

In  each  of  the  above  cases,  the  force  probe  and  sup¬ 
port  are  mounted  on  a  platform  suspended  by  rubber  supports 
which  in  turn  is  attached  to  a  linear-bearing  traversing  mechan- 
i sr  .  Model-orifice  alignment  is  carried  cut  using  optical 
techniques.  The  model  can  be  positioned  with  respect  to  the 
orifice  with  an  accuracy  of  0.2  mm. 
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3.3  Non-Ideal  Free  Jet  Flow  Field 


The  flow  properties  in  a  free  jet  do  not  always  cor¬ 
respond  to  the  "ideal"  picture  presented  above.  This  section 
discusses  the  practical  case  of  a  "non-ideal"  free  jet. 

There  are  two  classes  of  non-ideal  free  jet  flows 
which  need  to  be  considered:  those  which  are  due  to  the  manner 
in  which  the  free  jet  is  employed  in  this  particular  study 
(primarily  resulting  from  flow  moaulation),  and  those  which 
are  always  observed  when  free  jets  at  low  densities  are  employed 
for  aerodynamic  testing  purposes. 

Looking  at  the  first  class,  it  is  clear  that  modulating 
a  free  jet  with  a  rotating  chopper  disc  produces  a  periodic 
partial  blockage  of  the  flow  which  in  turn  causes  the  stagnation 
pressure  to  rise  and  fall.  Thus,  the  model  will  be  subjected 
to  a  tine-varying  force  when  the  orifice  is  uncovered  and 
covered.  This  problem  is  easily  circumvented  by  employing  a 
sufficiently  large  stagnation  chamber  such  that  the  response 
time  is  much  greater  than  the  disc’s  hal  b-per  i  od .  In  cold  flow 
studies,  a  50  liter  settling  chamber  is  employed  which  results 
in  pressure  fluctuations  of  a  few  tenths  of  one  percent.  The 
use  of  heated  flows  precludes  such  a  large  volume  because  of 
the  heating  requirements  that  would  be  imposed.  Pressure  fluctua¬ 
tions  in  the  one  liter  settling  chamber  used  for  hot  flows  were 
reduced  to  less  than  1*  by  positioning  the  chopper  disc  a  few 
millimeters  downstream  of  the  orifice  plane.  Flow  effects  felt 
on  the  free  jet  centerline  due  to  spillage  around  the  disc 
during  the  "flow-off"  period  were  investigated  with  impact 
probes  and  found  to  be  entirely  negligible. 

Flow  transients  must  also  be  expected  when  using  a 
chopper  disc.  Their  form  and  duration  have  been  studied  exper¬ 
imental!/  with  a  hot-wire  system.  When  the  flow  is  modulated, 
the  typical  response  of  a  nominal  one  micron  wire  is  that  shown 
on  Figure  7,  The  peaks  correspond  to  the  "closing"  of  the  ori- 
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fice  by  the  rotating  disk,  and  appear  to  be  quite  clean  and 
sharp.  The  valleys,  on  the  other  hand,  correspond  to  the 
"opening"  of  the  orifice  and  are  characterized  by  a  small 
irregularity  which  is  typically  2-3$  of  the  half-period.  It 
is  properly  termed  the  flow-establishment  time.  Note  that  the 
"uncovering-time"  of  the  orifice  is  approximately  1$  of  the 
half-period  (0.5  cm  dia,  orifice/50  cm  half-circumference  of 
the  rotating  disk).  Thus,  the  flow  transient  is,  es  might  be 
expected,  somewhat  in  excess  of  the  time  required  to  uncover 
the  orifice.  However,  the  ultimate  effect  of  this  transient 
on  the  force  measurements  is  diminshed  considerably  by  the  man¬ 
ner  in  which  the  lock-in  amplifier  treats  such  transients  (see 
Reference  2  for  details;  basically,  the  error  due  to  a  transient 
of  duration  At  is  not  Lt/i%  where  r  is  the  period  of  rotation, 
but  rather  (ht/x)’. 

Looking  at  the  second  class  of  problems,  it  is  seen 
that  when  using  the  results  of  the  method  of  characteristics 
(Ref,  3)  to  determine  the  flow  field  properties,  the  effective 
orifice  diameter,  D  ,  for  one  dimensional  flow  mu3t  be  known. 

It  is  nearly,  but  not  precisely,  equal  to  the  geometric  dia¬ 
meter.  For  each  orifice  employed,  the  effective  orifice 
diameter  was  calculated  as  a  function  of  the  Reynolds  number 
by  measuring  the  mass  flow  as  a  function  of  stagnation  pressure. 
For  all  sharp-edged  orifices  analyzed  to  date,  the  effective 
diameter  was  found  to  be  only  a  few  percent  smaller  than  the 
geometric  diameter  and,  more  importantly,  to  vary  by  no  more 
than  two  percent  over  the  entire  range  of  source  Reynolds  num¬ 
bers  employed  in  the  eventual  test  program.  For  purposes  of 

K 

computing  flow  properties,  D  was  assumed  constant. 

A  check  on  the  applicability  of  the  method-of-charac- 
teristics  solution  to  the  free  jet  was  rnRde  by  measuring  the 
impact  pressure  along  the  centerline  of  the  free  jet  as  a  func- 
tion  of  x/D  and  comparing  the  measured  values  with  the  pre¬ 
dicted  results.  Both  heated  and  unheated  flows  were  employed. 

The  raw  experimental  points  are  shown  in  Figure  6.  Corrections 
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for  shock  stand-off  distance  will  shi^t  all  the  data  points  to 
the  left  by  about  0 . U  x/D  .  Viscous  corrections  to  the  impact 
pressure  measurements  are  negligible  for  the  unheated  cases, 
but  they  will  shift  downwards  the  data  points  taken  under 
heated  conditions  by  as  much  as  6-8  %  (Ref.  *0  at  an  x/D*  of 
^12.  These  considerations,  combined  with  the  computed  measure- 
ment  accuracy  -  about  il%  at  an  :*/ D  of  5  ana  *8#  at  an  x/D  of 
12  -  allow  us  to  conclude  that  the  impact  pressure  surveys  are 
in  good  agreement  with  theory.  They  also  lend  credence  to  the 
fact  that  the  flow  is  not  disturbed  by  either  the  rotating  disk 
or  the  radiation  shields. 


A  final  non-ideal  free  jet  characteristic  is  the 
thickening  of  the  shock  waves  which  form  the  "barrel"  causing 
the  boundaries  to  encroach  on  the  useful  testing  region. 

Electron  beam  measurements  of  the  free-jet  density  profiles  on 
the  jet  centerline  have  been  carried  out  for  different  values  of 
the  stagnation  pressure.  The  results  are  shown  in  Figure  9.  The 
empirical  high-density  location  of  the  Mach  disk  (x^/D  )  which 
terminates  the  hypersonic  flow  field  is  calculated  from  the 
expression  (Ref.  3) 


•y  Po 


1/2 


D*  th 


(3-1) 


It  is  clearly  seen  that  reducing  po,  so  as  to  increase  X  ,  can 
only  be  continued  to  a  certain  level,  depending  on  the  axial 
location  of  the  model.  Below  this  level,  the  free  jet  density 
profile  no  longer  agrees  with  the  method  of  characteristics' 
predictions.  The  results  presented  herein  all  were  obtained  at 
values  of  po  such  that  the  measured  values  of  density  agree 
with  the  theoretical  predictions.  Therefore,  no  changes  in  flow 
properties  occur  due  to  the  encroachment  of  the  surrounding 
shock  structure  as  po  is  decreased  (so  as  to  increase  Kn  ), 


In  summary,  satisfactory  operating  conditions  for  the 
proper  use  of  a  free  jet  flow  field  have  been  confirmed  by  e 
variety  of  experimental  checks.  Within  the  boundaries  specified, 


the  theoretical  values  for  the  free  jet  flow  properties  may  be 
used  with  confidence  for  data  reduction  in  both  cold  and  hot 


flows,  requiring  only  a  knowledge  of  the  distance  from  the 
orifice  to  the  model  and  the  effective  diameter  of  the  orifice 
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k.  EXPERIMENTAL  RESULTS 

This  section  describes  the  test  procedure  and  the 
results  of  force  measurements  on  the  9°  half-angle  cones. 

A .  Test  Procedure 

The  modulated  force  probe  was  fitted  with  a  O.h  mm 
diameter  sting  which  supported  the  cone  model  in  the  manner 
shown  on  Figure  2.  The  orifice,  from  which  the  free  jet.  ema- 
nated,  had  an  effective  diameter,  D  ,  of  5 • ^  mm.  Based  on  the 
desired  Mach  number,  the  corresponding  axial  distance  from  the 
orifice  to  the  cone  was  calculated  according  to  the  free  jet 
theory  discussed  earlier  (Ref.  3).  The  sting  was  optically 
aligned  parallel  to  the  flow  centerline  as  determined  by  the 
plane  of  the  orifice. 

After  stabilizing  the  wind  tunnel  flow  and  the  elec¬ 
tronic  chain,  ‘the  stagnation  pressure  was  systematically  varied 
over  a  range  from  roughly  1  —  5 0  Torr  and,  at  each  value,  the 
force  was  measured.  Near  the  above  limits,  the  ambient  pressure 


was 

increased 

by  =25$  to  examine  its 

effect  on 

the 

readings;  if 

any 

effect  was 

noted,  the  reading  was 

discarded 

as 

having  been 

in  f  1 

uenced  by 

the  surrrounding  shock 

structure. 

When 

measurements  of  the  fo 

rce  acting 

on 

the  cone- 

sting  combination  were  completed,  the  cone  was  removed  from  the 
sting  and  suspended  by  a  fine  support  wire  just  in  front  of  the 
sting,  providing  the  proper  geometry  for  a  tare  measurement;  i.e. 
a  measurement  of  the  sting  drag  alone. 

A  calibration  was  performed  before  and  after  each 
wind  tunnel  series.  If  the  two  calibration  factors  differed  by 
more  than  2%,  the  data  were  discarded.  If  less,  the  average 
value  was  chosen  for  data  reduction  purposes. 
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B.  Results  in  Unheated  Flows 


The  force  measurements  at  a  given  pg  are  normalized 
by  the  appropriate  free  stream  properties,  computed  from  a 
knowledge  of  model  position,  x/D*,  and  pg,  to  form  the  desired 
dimensionless  grouping  -  the  cone  drag  coefficient  : 


CD  “  7 


DRAG 


2  pU»  AREF 


4—  1 


where  AREp  =  tD2/4.  The  distance  x,  is  taken  as  the  distance 
from  the  orifice  to  the  "half-area"  point  of  the  cone;  i.e., 
the  point  at  which  half  the  cone  area  lies  upstream  and  half 
downstream,  excluding  the  base  area. 


The  results  are  shown  in  Figure  10  for  different 
Mach  numbers  as  a  function  of  the  wall  Reynolds  number. 
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The  data  are  tabulated  in  Table  1  and  the  symbols 


are  identified  in  Table  2.  The  accuracy  of  the  data  was  computed 
based  on  the  following  considerations  : 
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A  best  estimate  of  the  random  error  is  tabulated  below 


M  *6 


^8 


*10 


Force 

Po 

x 

D* 

A 

REF 


±1.0# 

±1.0# 

±1 . 0#-±0 . 05# 

±0.2#- 

±0.5# 

±0.4# 

±0.5# 

±0.5# 

±1.0# 

±1.0# 

±1.0# 

04#  ±0.2#-±0.04# 

± 0.2 # 

±0.5# 

±1.0# 


R .  M .  S . 

Random  ±2.2#-±2.0#  ±1.9#  ±1.8# 

Error 


Thus  a  random  error  of  approximately  ±2#  in  can  be 
expected.  Naturally,  the  presence  of  systematic  errors  is  also 
possible.  In  particular,  the  Mach  6  data  appear  to  be  low  by 
roughly  15#i  this  point  is  discussed  in  a  later  section. 


The  data  presented  in  Fig.  10,  recorded  in  an  unheated 
flow,  include  the  effects  of  half-angle,  0c;  type  of  gas,  y;  and 
Mach  number,  M^.  A  qualitative  test  on  the  effect  of  bluntness 
vas  carried  out  at  Mach  10.4  with  two  cones  having  bluntness 
ratios  R„/Rr>  0.04  and  0.30.  When  plotted  versus  Re  for 

example,  the  blunt  cone  exhibited  slightly  lower  (  =  5#)  values 
of  C  ^  than  did  the  sharp  cone. 

C.  Results  in  Heated  Flows 


As  tests  in  a  heated  flow  were  expected  to  raise  the 
overall  temperature  level  of  the  force  probe  by  a  not-insigni- 
ficant  amount,  calculations  were  first  performed  to  estimate 


The  larger  value  applies  to  the  lowest  value  of  pg  (and  thus 
the  lowest  Re  ) ,  the  smaller  value  applies  to  the  largest 
value  of  p0  (and  thus  the  highest  Re  ). 

V  y  X, 
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the  equilibrium  temperature  of  the  cone  model  under  a  variety 
of  flow  conditions.  Factors  included  in  the  calculation  were  : 
radiation  heat  transfer  from  the  interior,  exterior  and  surround¬ 
ing  shields  of  the  settling  chamber;  convective  heat  transfer; 
and  conduction  heat  transfer  through  the  connecting  sting  and 
support  assembly.  Detailed  results  can  be  found  elsewhere 
(Ref.  5).  Of  essential  interest  was  the  fact  that  model  temper¬ 
atures  between  350~550°K  might  be  expected  with  the  existing 
probe  design  for  the  anticipated  range  of  flow  conditions  and 
a  stagnation  temperature  of  1200°K  -  the  maximum  design  temper¬ 
ature  for  the  heated  stagnation  chamber  discussed  earlier. 
Considering  the  specific  geometry  of  the  force  probe,  it  was 
further  concluded  that  the  transducer  temperature  would  be 
about  330°K  50°C).  This  is  well  below  the  Curie  temperature 

of  the  PZT-5  ferroelectric  employed,  but  it  is  known  that  the 
properties  of  these  crystals  are  somewhat  dependent  on  temper¬ 
ature  over  such  a  range. 

A  calibration  of  the  transducer  was  therefore  performed 
in  a  temperature-controlled  oven  at  temperatures  up  to  360°K. 
Small  1055)  changes  in  the  calibration  factor  of  the  ferro¬ 
electric  device  were  noted,  but  the  performance  was  otherwise 
quite  stable  and  repeatable. 

When  tests  in  a  heated  flow  in  the  wind  tunnel  were 
initiated,  however,  completely  different  results  were  obtained 
characterized  by  marked  instabilities  and  force  values  which 
were  far  removed  {'\5%-25%)  from  the  expected  values.  Measurements 
of  the  model  temperature  confirmed  the  earlier  predicted  values 
and  hence  the  transducer  temperature  vas  certainly  well  within 
the  limits  reached  in  the  oven  during  calibration  (i.e.,  <360°K). 
The  only  essential  difference  between  the  oven  calibration  tests 
and  the  wind  tunnel  experiments  was  that  in  the  latter  a  temper¬ 
ature  gradient  existed  along  the  length  of  the  transducer  (a  AT 
of  approximately  10°C),  whereas  the  former  were  performed  at  a 
uniform  element  temperature. 
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Since  it  was  desirable  to  further  reduce  the  model 
temperature  to  provide  an  increased  ratio  of  To/Tw»  and  as  such 
a  change  would  also  decrease  both  AT  across  the  transducer  as 
well  as  the  average  temperature  of  the  transducer,  the  probe 
was  modified  to  provide  an  alternate  conduction  path  for  the 
heat  which  would  reduce  the  model  and  transducer  temperature 
levels  at  the  expense  of  only  a  small  loss  in  sensitivity.  The 
result  was  unsuccessful  in  the  sense  that  the  same  unstable 
operation  was  observed. 

Finally,  a  special  test  set-up  was  arranged  to  quanti¬ 
tatively  examine  the  effect.  In  the  calibration  mode,  the 
supports  fore  and  aft  of  the  transducer  (see  Fig.  3)  were  fitted 
with  thermocouples  and  a  blackened  foil  was  cemented  with 
conducting  glue  to  the  sting.  The  light  from  a  150  watt  reflector 
lamp  was  focused  onto  the  black  foil.  Thus,  sufficient  heat  was 
added  forward  of  the  transducer  to  produce  a  AT  across  the 
transducer  which  could  reach  30°C.  At  the  same  time,  a  force 
could  be  applied  to  the  element  and  thus  the  effect  of  a  well- 
defined  temperature  gradient  on  the  output  signal  could  be 
accurately  determined.  An  oscilloscope  provided  a  direct  record¬ 
ing  of  the  signal  at  the  input  to  the  lock-in  amplifier.  The 
results  were  striking.  A  heat  pulse,  producing  an  initial  rate 
of  change  of  AT  on  the  order  of  0.1°C/sec,  caused  the  transducer 
output  to  surge  in  such  a  way  as  to  saturate  the  charge  amplifier 
and  momentarily  cut  to  zero  the  input  signal  to  the  lock-in 
amplifier.  Some  tens  of  seconds  were  required  for  the  system  to 
return  to  normal  operation.  Various  tests  eventually  proved 
conclusively  that  the  responsible  cause  was  the  change  in  element 
temperature  with  time.  The  mechanism  is  probably  a  little-under- 
stood  (from  the  quantitative  point  of  view)  phenomena  called 
the  pyro-electric  effect.  It  should  be  emphasized  that  the  pyro¬ 
electric  effect  can  easily  produce  volts,  whereas  the  largest 
output  signals  observed  in  force  measurements  are  typically 
millivolts.  Thus,  complete  saturation  is  quite  understandable. 
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When  a  ti me- i ndependent  AT  of  2°C  was  established  across 

the  element,  the  output  signal  was  quite  satisfactory.  At  a  AT 
of  k°C,  some  unstable  operation  was  noted;  enough  to  remove  the 
element  from  the  realm  of  a  ±1#  instrument.  At  larger  AT's, 
operation  was  impossible  due  to  large  drifts  and  instabilities. 

Since  the  applied  force  during  all  of  the  above  tests 
was  relatively  large  compared  to  the  values  recorded  under  near- 
free  molecule  flow  conditions,  it  became  clear  that  even  a 
AT  =  2°C  could  not  be  tolerated. 

In  summary,  the  heated  flow  tests  could  not  be  carried 
out  with  any  reasonable  degree  of  accuracy  (<  2055)  because  of 
a  highly  unstable  output  signal  produced  by  a  temperature 
gradient  of  only  a  few  degrees  Celsius  or  more,  a  gradient  so 
small  that  any  modifications  to  the  probe  design  for  the  purpose 
of  providing  increased  cooling  would  have  to  be  so  important  as 
to  produce  unacceptable  side-effects;  i.e.,  a  highly  reduced 
sensitivity  due  to  the  mechanical  constraint  of  the  heat-conduction 


member . 
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5.  DISCUSSION  OF  RESULTS 

An  examination  of  Figure  10  discloses  a  number  of 
points  which  are  commented  on  here. 

1 .  The  Mach  8 . 1  data  <x,+  ) ,  while  exhibiting  little  scatter 
within  each  set,  show  a  systematic  difference  from  each  other 
which  varies  from  6  to  8%.  A  similar  result  is  noted  for  the 
Mach  10.  U  data  (□ ,  a  ) •  In  principle,  the  data  for  each  Mach 
number  should  be  averaged  resulting  in  the  conclusion  that 
for  these  two  sets  of  tests  the  reproducibility  is  roughly 

i  3  to  1*#.  As  the  tare  correction  typically  amounted  to 
30  -  Uo$3  of  the  total  drag,  some  discrepancy  may  be  due  to  the 
inevitable  substraction  errors  that  occur  in  a  sting-mounted 
arrangement.  A  second  source  of  error  is  in  the  location  of 
the  model  relative  to  the  orifice.  Because  of  the  expanding 

nature  of  the  flow,  the  Mach  number,  and  hence  normalized 

* 

dynamic  pressure  and  changes  with  yj D.  As  an  example,  an 
error  in  x  of  V2  mm  will  result  in  a  2  -  3%  error  in  the 
calculated  value  of  for  the  Mach  8  case. 

2.  While  the  bulk  of  the  data  falls  on  a  straight  line 

emanating  from  the  C/C  =  1  point,  a  tendency  is  noted 

D  D  f .  m . 

for  the  data  to  rise  above  this  line  for  Rew>£<  2  or  3.  This 
fact  has  been  noted  previously  (Ref. 1,  p.33),  but  whether  the 
effect  is  real  or  simply  represents  experimental  difficulties 
in  the  regime  of  very  low  force  levels  is  not  clear. 

3.  The  data  have  been  compared  with  the  predictions  of  Willis' 
Knudsen  iteration  analysis  (Ref.l)  which  provides  a  first  order 
correction  to  the  free  molecule  drag  coefficient  and  is  based 
on  the  modified  Krook  model  for  the  collision  term  in  the 
Boltzmann  equation.  The  expansion  parameter  in  Willis'  analysis, 
essentially  a  Knudsen  number,  can  be  defined  in  different  ways 
depending  on  the  test  conditions.  Briefly,  Willis  predicts 
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where  A  is  a  function  of  cone  angle  0  ,  and  wall  speed  ratio, 

c 

5  ,  and  a,  is  taken  to  be 
w  oc 


be 
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where  Kn ( D) 


*0p  o 

D  p_ 


When  these  auantities  are  rewritten  in  terms  of  Re  the 

w ,  l 

straight  lines  on  Figure  10  result. 


It  is  noted  that  the  agreement  between  theory  and 

experiment  is  very  good  for  the  0  =  13*5°  case  and  reasonably 

c 

good  for  the  0  =9°  case  (with  the  exception  of  the  Mach  6 

c 

case  discussed  below),  particularly  when  slopes  rather  than 
absolute  values  are  considered. 


.  The  Mach  6  data  are  appreciably  below  (=  15$)  the  other 

data.  While  some  Mach  number  effect  is  anticipated,  such  a 

discrepancy  is  too  large  to  be  expected.  Noting  that  a  straight 

line  through  the  Mach  6  data  is  nearly  parallel  to  the  theoretical 

line  (the  solid  line  for  G  =  9°),  it  is  possible  that  a  systematic 

c 

error  occurred.  Extrapolating  the  data  to  Rew  0  yields  a 

C ,/Cr,  =  0.86  instead  of  unity.  If  the  values  of  the  ordinates 

^ •  rn  ♦ 

are  raised  by  a  corresponding  amount,  then  the  "corrected" 
data  for  the  Mach  6  case  will  (naturally)  fall  on  or  around  the 
solid  line  and  at  the  same  time  will  agree  with  the  other  data. 
Clearly,  no  firm  theoretical  or  experimental  justification  exists 
for  this  type  of  data  manipulation.  Such  a  large  systematic 
error  could,  however,  arise  due 'to  an  error  in  the  measurement 
of  x,  the  distance  from  the  orifice  to  the  model.  In  this  case, 
an  error  of  1.5  mm  would  be  necessary  to  explain  the  observed 
discrepancy.  A  more  likely  explanation  is  related  to  the 
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fundamental  question  regarding  the  selection  of  x  itself. 

Is  it  the  distance  to  the  apex,  the  mid  point,  the  "half-area" 

point  (used  herein),  or  the  base  ?  As  the  models  have  a 

nominal  length  of  5  mm,  it  is  readily  seen  that  the  normalized 

dynamic  pressure  (and  hence  the  calculated  C^)  will  vary  by 

=  from  apex  to  base  for  the  Mach  6  case  because  the  distance 

^  C  $  x 

from  orifice  to  apex  is  21.  «  mm  and  - — -  a,  2  .  (The  variation 

C  D  X 

of  the  normalized  dynamic  pressure  for  the  Mach  8  to  10  cases  is, 
of  course,  much  less  because  x  is  larger;  the  corresponding 
percentages  would  be  =  25$  and  -  15$  respectively.)  As  it  ir 
quite  difficult  to  accurately  machine  smooth  slender  cone  moi  -Is 
with  a  length  of  less  than  5  mm,  and  as  the  absolute  size  of 
the  VKI  free  jet  is  quite  large  compared  with  others  currently 
in  use  for  similar  studies,  the  wide  variation  in  calculated 
values  of  which  are  possible  depending  on  the  choice  of  x 
should  be  kept  in  mind  wnen  analyzing  free  jet  results.  In 
summary,  although  an  evaluation  of  flow  properties  at  the  "half¬ 
area"  point  appears  logical,  it  cannot  be  said  to  be  correct 
in  an  absolute  sense. 


19 


6.  A  CORRELATION  PARAMETER  FOR  CONE  DRAT. 


It  is  clearly  of  interest  to  develop  a  parameter 

vhich  would  correlate  cone  drag  coefficients  over  a  wide  range 

of  conditions.  In  the  past,  experimentalists  in  the  low  density 

field  generally  employed  the  free  stream  Knudsen  number  Kn^, 

while  members  of  the  continuum  school  used  the  free  stream 

Reynolds  number,  Re  .  Ignoring  the  effect  of  relative  velocity 

00 

on  collision  frequency,  the  two  quantities  are  related  in  a 
simple  way  as  first  shown  by  von  Karmen;  Kn  a,  M  /Re  where 

03  00  00 

the  coefficient  is  of  order  unity.  Thus  the  two  parameters  are 
not  at  all  independent,  but  on  the  other  hand  the  success  of 
either  one  is  limited  as  can  be  seen  in  Figure  11,  for  example, 
where  a  wide  variety  of  data  has  been  plotted  as  a  function  of 
Kn  *.  It  should  not  be  surprising  that  Kn  (or  Re  )  is 

CO  ^  Jj  00  00 

unsuccessful  in  the  low  density  regime  because  it  says  nothing 
about  the  effect  of  the  body  or  wall  temperature;  an  effect 
which  is  clearly  important  in  the  free  molecule  regime  and  thus 
one  which  must  persist  at  least  to  some  extent  as  Kn^  decreases. 

However,  there  is  general  agreement  that  a  Knudsen 

number  based  on  the  mean  free  path  of  a  molecule  leaving  the 

body  in  question  and  colliding  with  an  incident  or  free  stream 

molecule  (defined  here  as  X  . )  must  characterize  the  near-free 

n 

molecule  regime.  It  is  a  Knudsen  number  vhich  accounts  for  both 

the  free  stream  properties  (M  ,  T^)  and  the  properties  of  the  wall 

(T  )  and  may,  in  addition,  account  for  the  temperature  (or 
w 

velocity)  dependence  of  the  cross  section  characterizing  this 
particular  type  of  collision.  The  latter  point  seems  important 
to  this  author  when  comparing  results  in  which  the  relative 


W’U  U  n  n  o  2  ft  (  T  ) D 
00  00 

where  fi(T  )  may  be  calculated  from  tables  presented 

CO 

in  R e f . 6 . 
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velocities  characterizing  collisions  in  the  various  wind  tunnel 
experiments  on  record  vary  by  nearly  one  order  of  magnitude. 
Such  a  Knudsen  number  has  been  successful  in  correlating  a 
wide  variety  of  sphere  drag  data  (Ref. 7). 

The  purpose  of  this  section,  then,  is  to  define  a 
Knudsen  number  based  on  a  reflected-incident  collision  process 
which  includes  a  realistic  temperature  (or  energy)  dependent 
cross  section  and  then  to  use  this  so-called  "effective" 

Knudsen  number  to  correlate  a  variety  of  cone  drag  results 
and  assess  the  consequences.  This  section  essentially  repeats 
material  previously  presented  in  Ref. 7- 

The  definition  of  the  "effective"  Knudsen  number  and 
the  definition  of  the  mean  free  path  associated  with  the 
collision  process  is  : 


Kn 


ri  ,D 


X  . 
r  1 

D 


D 


where  | g  . |  =  | ( £ .  -  X  )  |  • 

i  -rl  i  1  *r  ' 


6-  I 


£ince  a  near-free  molecule  analysis  indicates  that 
the  first  order  correction  to  the  free  molecule  drag  coefficient 
(see,  e.g.,  Ref. 8)  goes  as  : 


' D „  CD  %  Kn  . 

f .  m .  r  l 


therefore 


C 

D 


C 


D 


1 


f .  m. 


-v  1 


f .  m. 


Kn  . 
ri 
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This  point  is  emphasized  because  it  is  common  (and  indeed 
essential)  to  normalize  near  free  molecule  drag  by  the  free 
molecule  drag  coefficient  and  thus  the  relevant  correlation 
parameter  is  more  correctly  C  Kn  .  rather  than  just 

Df.n.  ri 

Kn  •  .  (Since  C_  seldom  exceeds  3,  and  then  only  at  los 
?f,m. 

supersonic  conditions,  and  cannot  be  less  than  2,  the  inclusion 
of  this  parameter  as  a  multiplying  factor  does  not  exert  an 
important  influence  on  the  results). 

It  remains  to  put  equations  6-1  into  a  form  which 
will  allow  the  insertion  of  experimentally  measurable  quantities. 
Assume  that  we  can  write  : 


Choosing  the  x-direction  parallel  to  the  bulk  velocity  U^; 
putting  the  cone  at  the  origin  of  the  coordinate  system  so 
that  molecules  reflected  from  it  have  no  bulk  velocity,  but 
only  their  random  thermal  velocity,  so  that  X  =  X  ;  and 

3F-  r  r 

using  the  fact  tnat  C  =  0;  we  have 

00 


(f.  -  I  )2  =  U  2  +  2U  C  +  C  2  +  C  2 

l  x*  oo  ®  x  oo  X* 

X 

Ignoring  the  complicating  effects  of  geometry  and  simply 
picturing  all  molecules  as  emanating  from  a  plane  perpendicular 
to  the  free  stream  velocity. 


(It  is  clear  that  this  assumption  is  less  adequate  for  cones 
than  for  blunt  bodies;  however  the  quality  of  experimental 
data  to  date  probably  does  not  justify  additional  rigor.) 
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Recalling  that  C2  =  3RT,  the  main  term  of  interest  in  6-1 
becomes 


|  5 


|<VU  /  <s, 


2  ♦ 
3 


“  1  +  ~ 

S2  J  3 


2  S 


2*  +  1) 
/it 


6-14 


(!!ote  that  for  S  >>  1  ,  S  >>  1;  6-h  is  nearly  S  ^). 

oo  w  w 


We  define 


F  ( S  ,S  )  =  S  2  [  |  +  - 

w  *  v  [  3 


S2 


2  Sw 

+  ~  —  +  1 


3  /, 
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Therefore , 

Kr. 


X  . 
ri 


rl,D  D  n  Z(r  .  )  D  /F  (  ^  ,S_) 

m  it  1  rr. 


®  r  i 


w  ® 


6-6 


The  velocity  dependent  cross  section  E(g^.)  must 
now  be  put  into  a  sensible  form.  We  compute  an  equivalent 
temperature  characterizing  the  collision  process  by  computing 
the  temperature  of  an  equilibrium  gas  in  which  the  mean  square 
relative  velocity  between  two  molecules  is  taken  equal  to  tne 
mean  square  relative  velocity  of  a  free  stream  molecule  with 
respect  to  a  molecule  emitted  from  the  surface;  i.e.  in  a 
general  equilibrium  case  : 


g2  =  C,2  +  C22  =  2 ( 3  R  T  ) 

1  *■  e  q 

In  the  situation  described  here. 


(t.  -  £  )2  =  U2  +  2U  c 
i  r  oo  oo 


r 


x 


t-  c2  +  c2 , 

oo  r 
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which  is  the  mean  square  relative  velocity  for  the  collision 

geometry  considered.  Equating  the  two  allows  one  to  define  an 

equivalent  temperature  for  the  collision  process  which  is 

defined  in  turn  as  T  .  and  is  given  by  : 

r  1 


2 


T  F  ( S  ,S  ) 
w  v  * 


6-7 


where  F(S  ,S  )  was  defi.ed  in  6-6. 
w  « 

Note  that,  as  expected,  when  >>  1  and  »  1  (cold  wall 
hypersonic  case )  , 


T  . 
r  l 


1  . 1 . 

3  y  -  i 


To  *  T0 


Thus,  the  effective  Knudsen  number  characterizing 
collisions  between  free  stream  molecules  and  molecules 
emanating  from  a  plane  at  a  temperature  T  placed  perpendicular 

V 

to  the  flow  direction  may  be  computed  from  equations  6-6  and 
6-5  where  the  cross-section  is  computed  at  the  temperature  T  ^ 
given  by  equation  6-7- 


It  is  instructive  to  note  the  relationship  between 

Kn  .  and  parameters  employed  by  other  investigators.  A  few 
r  l 

such  comparisons  follow. 


I.  Consider  the  case  :  F  >>  1 ,  F  >  S  (the  hypersonic  cold- 

w  «  w 

wall  case).  From  equation  6-5,  F(Rw,  )  -*■  S2  and  T  ^  T0 . 
Then 


Kn  . 


Kn  Z(T  ) 

CO  00 


ri,D  n(oI(To)DSw  Sy  Z(T0) 


Kn 


If  E  =  constant  (hard  sphere  case),  then  Kn  .  _  %  - — , 

r  l  ,  v  b 

, 

a  parameter  used  m  Ref . 1 ,  for  example. 

If  E  ^  T~ 1 / 2  (implying  y  ^  T;  thus,  a  "Krook"  gas),  then 


Kn  S 

CO  cc 

Kn  .  _  n,  — - - 

r  i  ,  D  S 
’  w 


Since  y  a,  Tl/2/'E(?),  then  using  the  fact  that  Kn  n,  M/Re, 


the  correlation  employed 


it  is  seen  that  Kn 
in  Ref.  9* 


ri  ,D 


1 

Re0 
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II.  Next,  consider  the  case  :  S  a,  1  (the  hypersonic  hot  wall 

case).  From  equation  6-5,  F ( S  ,  S  )  +  1  as  long  as  S  >  S 

w  “  ®  -  w 

and  T  .  -*■  T  =  T0  when  S  >>  1.  Then 
n  w  u  ® 


Kn 


ri  ,D 


1 

n  Z  (T  )D 

ot>  V 


Kn 

00 


z  (T  ) 

GO 


If  £  =  constant,  then  Kn  .  „  e  Kn  . 

ri  ,v  ® 

At-1 

If  I  ^  T  J/2,  then  Kn  .  _  Kn  /  ~  ;  but  since  T  -  Tn 

r  l  ,  D  »  /  T  ’  w  0 

when  S  ^  1  if  S  >>  1  ,  it  follows  that  Kn  .  _  Kn  M  .  Thus, 
w  “  ri  ,D  a>  m 

it  is  readily  seen  that  the  effective  correlation  parameter  for 

this  problem,  restricting  our  attention  to  the  hypersonic  case, 

may  be  expressed  as  the  product  of  Kn^  and  a  quantity  which 

will  vary  from  1/S^  to  S^/R^  to  1  to  Sm  depending  on  the 

temperature  dependence  of  the  collision  cross  section  and  the 

relative  magnitude  of  S  compared  to  unity. 

w 

It  is  seen  that  the  effective  Knudsen  number  defined 
in  6-1  reduces  to  the  conventional  forms  of  Kn  (or  Re)  under 
certain  limiting  conditions.  Therefore,  it  is  clear  that  a 
realistic  temperature-dependent  cross  section  must  be  incorporated 
ir  any  such  correlation  paramo  if  a  wide  variety  of  data 
representing  many  distinctly  different  flow  conditions  are  to 
be  adequately  correlated. 

A  limited  amount  of  slender  sharp  cone  drag  data 

has  been  plotted  in  terms  of  Kn  .  „  Cn„  on  Figure  12.  The 

ri ,D  wf . m. 

symbols  are  defined  in  Table  2.  The  wall  speed  ratios  vary  from 

1.5  to  10,  y  from  1.30  (C07)  to  1.67(A),  M  from  3*7  to  2b , 

00 

and  0  from  6°  to  15°- 
c 
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A  word  of  caution  must  be  expressed  (and  possibly 
an  apology  as  well)  regarding  the  translation  of  data, 
originally  expressed  in  terms  of  Kn^,  Re^ ,  etc.,  into  Kn^. 

In  some  cases,  a  bit  of  guessing  as  to  the  values  of  unspecified 
flow  or  model  conditions  was  required  to  effect  the  transformation. 
In  other  ca^es,  certain  flow  conditions  were  recomputed;  e.g», 
Geigers  values  of  Rew  in  Ref. 10  were  recalculated  because 
the  Sutherland  viscosity  law  used  therein  is  well  known  to 
be  incorrect  at  values  of  T  below  50°K.  The  Chapman-Ensk og 

CO 

viscosity  with  temperature-dependent  cross  sections  (Ref. 6) 
was  employed  in  all  data  manipulation  for  consistancy.  In 
addition,  the  data  of  Kussoy  et  al  (Ref. 12)  was  actually  taken 
from  the  graphs  of  Ref.  1. 

The  results  in  Figure  12  are  not  as  strikingly  clear 

as  in  Ref. 7  for  sphere  urag  where  much  more  data  were  available. 

In  that  case,  it  was  found  that  data  for  different  y's  and 

M  's  fell  on  a  single  curve  for  a  given  value  of  S  and  that 
»  w 

as  S  increased,  C^/Cr,^  increased  for  a  fixed  value  of  Kn  .  n. 
w  D  f.m.  r  l  ,  u 

Here  the  same  trend  can  be  inferred  only  if  Geiger's  data  are 
ignored.  It  should  be  mentioned  that  Willis'  theoretical 
predictions  (Ref. l)  lie  above  Geiger's  results,  below  the  results 
of  Jones  and  Miller  (Ref. 11),  and  generally  somewhat  below 
the  results  of  Kussoy,  et  al.  (Ref.  12).  Willis'  theory  also 
predicts  the  same  variation  of  C^/Cn^  with  S  as  suggested 
here . 


A  close  inspection  of  Figure  12  illustrates  that  the 

quantity  Kn  .  _  is  a  useful  correlation  parameter, 

r  l  ,  D  .  D  f .  m . 

This  assertion  is  borne  out  when  three  sets  of  data  for  nearly 

the  same  are  plotted  in  Figure  13.  Note  that  although  y 

varies  from  1.3  to  1.6'7  and  M  varies  from  3-7  to  11.5,  a 

00 

satisfactory  correlation  is  achieved.  It  might  even  be 
argued  that  the  CO2  data  of  Sims  (Ref.  13)  would  be  expected  to 
lie  somewhat  above  the  other  data  because  the  cone  half-angle 
in  Sims'  case  is  somewhat  smaller  (10°)  than  the  other  cases 
(13.50)  _  again  a  trend  predicted  by  Willis'  theory. 
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7.  SUMMARY 

Drag  coefficients  for  sharp  9°  cones  have  been 
experimentally  determined  in  the  hypersonic  hot-wall  regime 
at  nominal  Mach  numbers  of  6,8  and  10.  The  results  are  in 
satisfactory  agreement  with  the  data  of  other  investigators 
and  with  Willis'  theory  for  the  near-free  molecule  regime. 

Cones  with  a  bluntness  ratio  of  0.3  exhibit  slightly  smaller 
values  of  at  otherwise  identical  conditions. 

The  correlation  parameter  Kn^.  ,  successfully  employed 
foi  blunt  bodies, appears  to  correlate  data  in  the  transitional 
and  near-free  molecule  regimes  for  varying  Mach  numbers  and  specific 
heat  ratios  if  the  wall  speed  ratio  is  constant. 

More  precise  data  on  slender  cones  are  required  in 
the  cold-wall  regime  before  a  more  definitive  statement  can  be 
made  regarding  the  adequacy  and  limits  of  applicability  of 
this  parameter. 
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TABLE  1 


Cone  Drag  Coefficients  -  VKI  Data 


A.  Half  Angle  -  9°  (sharp),  T0  =  T  =  293°K 
Gas-Air,  w 


Symbol 

Re  « 

w ,  *• 

Kn  „ 

00 ,  D 

c  /c 
d'  d 

* 

2.97 

1 .20 

0.820 

4.22 

0.850 

0.795 

M„  = 

6 . 0 

5.24 

0.684 

0.765 

Sw  = 

1.75 

6.60 

0.542 

0.745 

7-78 

0.460 

0.735 

11.1 

0.322 

0.685 

11.8 

0.304 

0.675 

13.5 

0.265 

0.665 

17.4 

0.206 

0.615 

21.5 

0.167 

0.580 

25.3 

0. 1 4 1 

0.390 

31.4 

0.114 

0.365 

41.6 

0.0860 

0.340 

52.2 

0.0685 

0.320 

y 

3.64 

0.725 

0.975 

A 

L .  94 

0.564 

0.920 

= 

8.1 

6.05 

0.460 

0.900 

Sw  = 

1.80 

7.15 

0.388 

0.875 

8.85 

0.314 

0.850 

11.6 

0.239 

0.790 

14 .7 

0. 189 

0.755 

JU 

3.03 

0.916 

1 .04 

I 

3.78 

0.736 

1 .02 

Moo  = 

8.1 

4.87 

0.571 

0.955 

Sw  = 

1.80 

5.84 

0.476 

0.950 

8.82 

0.315 

0.920 

11.6 

0.240 

0.850 

14.7 

0.189 

0.820 

□ 

1.78 

1.28 

1 .02 

2.07 

1.10 

0.990 

M«  = 

10.4 

3.42 

O.665 

0.910 

Sw  = 

1.83 

4.31 

0.527 

0.925 

B 

1.84 

1.23 

1.01 

2.36 

0.963 

0.925 

Moo  = 

10.4 

3.25 

0.698 

0.830 

Sw  “ 

1.83 

4.25 

0.535 

0.860 

2.34 


2.22 


2.22 


2.19 


2.19 


30  - 


B.  Half  Angle  -  13-5°  (sharp),  T0  =  Tw 
Symbol  Hey  £  Kn®,D  °D/CD 


=  293°K 


2.39 

0.808 

0.965 

2.42 

0.798 

0.975 

8.1 

3.12 

0.619 

0.940 

1.80 

3-77 

0.511 

0.915 

gen 

5.25 

0.368 

0.875 

5.26 

0.367 

0.890 

7.65 

0.252 

0 . 805 

11.6 

O.166 

0.720 

18.5 

0.104 

0.655 

3.92 

0.231 

0.960 

5.13 

0.177 

0.900 

11.5 

5.85 

0.155 

0.875 

1.56 

6 . 60 

0.137 

0.850 

1 

6.71 

0.135 

0.835 

7.31 

0.124 

0.825 

8.13 

0.111 

0.790 

9.62 

0.0943 

0.760 

13.0 

0.0695 

0.700 

18.9 

0.0480 

0.605 

18.9 

0.0480 

0.595 

34  •  1 

0.0266 

0.470 

35.8 

0.0253 

0.450 

51-9 

0.0175 

0.395 

2.25 


2.27 


TABLE  2 


Experimental  Conditions  for  Data  Presented 
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FIGURE  6  -SCHEMATIC  OF  INCONEL  HEATER 
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Fig. 8  IMPACT  PRESSURE  VARIATION  IN  FREE  JET 
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